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Alberto Lusiani 



1. Introduction 

In recent years, the BaBj^ and Belle experiments have contributed results on tau lepton physics, 
which improved the experimental picture of lepton universality and of lepton flavor violation 
searches. 

Both experiments rely on "B-factories" operating at a centre-of-mass energy of 10.58 GeV 
on the T{As) peak, just above the threshold for producing B-mesons. BABMi operates at the PEP- 
II complex at SLAC, which collides 9 GeV electrons against 3.1 GeV positrons, and has recorded 
about 420fb^^ of data by May 2007. Belle operates at the KEKB B-factory in Japan, which collides 
8 GeV electrons against 3.5 GeV positrons, and has recorded about VlOfb^' of data by May 2007. 

The BaBml ^ and Belle [§] detectors share several similarities and both include a silicon mi- 
crovertex detector, a drift chamber, a 1 .5 T solenoidal superconducting magnet, an electromagnetic 
calorimeter based on Cesium Iodide crystals, and a segmented muon detector in the magnet re- 
turn yoke. The two experiments differ in the particle identification strategy: Belle uses an aerogel 
threshold Cherenkov detector together with time-of-flight and tracker dE/dx, whereas BaBj^ relies 
on a ring-imaging Cerenkov detector supplemented by the dE/dx in the trackers. 

With a total now exceeding 1.1 ab^' of integrated luminosity and a e+e^ t+t^ cross- 
section at 10.58 GeV of 0.919 nb [Q], B-factories recorded in excess of 10^ tau pairs, which allow 
for improving statistics-hmited results, hke in particular searches for tau lepton flavor violating 
decays. 



2. Lepton universality tests 

The Standard Model (SM) predicts that all lepton doublets have identical couplings to the W 
boson. Ratios of measured decay widths of leptonic or semileptonic decays which only differ in 
the lepton flavour test whether the W interaction is universal for the three lepton flavours. The 
present data are consistent with the universality of the leptonic charged-current couplings to the 
0.2% level [Q]. B-factories have produced results on the tau mass, the tau lifetime, and could in 
principle measure tau branching fractions that improve the experimental knowledge on the less 
precisely known factors in the following expressions for leptonic coupling ratios, which the SM 
predicts to be consistent with unity: 

— = -^^{r -^e VeVr) — . , ^ (2.1) 



..... ; ' "2 7 (2-2) 



(2.3) 



In the above expressions, Y represents partial widths, lepton lifetimes, lepton masses, ^ 
branching fractions; f{x) = 1 — 8x + 8x^ —x^ — \2x\nx are phase space factors [||], /^-^ « 1 corre- 
spond to electro-weak radiative corrections [Bl]. 



2 



Lepton Universality and Lepton Flavour Violation tests at the B-factories 



Alberto Lusiani 



CM 
U 




500 - 



1.72 1.73 1.74 1.75 1.76 1.77 1.78 1.79 1.8 

M„,„ (GeV/c^) 

Figure 1: Pseudomass distribution forMmin for the T ^ 3;rv candidates. The tau mass is a parameter of the 
fit shown with a sohd fine. 



2.1 Tau mass 



The mass of the tau lepton appears at the fifth power in the coupling ratio expressions 2. 1 and 
2.2, contributing a relative uncertainty of 0.08% [|6|]. 



BELLE measured the tau mass using a pseudomass technique that was first employed by the 
ARGUS collaboration [Q]. Tau pairs candidates are selected where one tau decays into a single 
prong (electron or muon), and the other one decays into 3-prongs, all of which are pions, with no 
additional tt". Then a quantity M^m is computed: 



Mn,i„ = ^^Ml + 2(£beam -Ex){Ex-Px), (2.4) 

which is less than or equal to the tau lepton mass. Mx, Ex and Px are the invariant mass, energy and 
absolute value of the momentum, respectively, of the hadronic system in the center-of-mass (cm.) 
frame, and ^beam is the energy of the electron (or positron) in this frame. The distribution of Mmin 
extends up to and has a sharp edge at M^, smeared by detector resolution and initial and final state 
radiation. The edge position obtained from a fit to the M^am distribution is relatively insensitive to 
the background and is a precise estimator of the tau mass. The offset between the edge and the tau 
mass has to be estimated with a Monte Carlo. The measurement (on 414fb^^ of data) reaches a 
precision close to the present BES-dominated world average, = 1776.991q 26 

mi: = 1776.61 ±0.13 ±0.35MeV [§. (2.5) 

The limiting systematic contributions to this result come from the understanding of the mo- 
mentum scale and the reliability of the simulation in estimating the offset between the tau mass 
and the edge position. The uncertainty on the momentum scale is assessed by comparing the re- 
constructed B meson masses to their world average, which is know up to about 1.5MeV [Q], and 
noting that the di-muon invariant mass peak matches the center of mass event energy determined by 
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the beam energies withing 3 MeV. Apparently, measurement conducted at threshold with machines 
where the beam energy can be calibrated through resonant depolarization, suffer from smaller sys- 
tematic uncertainties, as the recent KEDR result shows (I776.8O+U 23 ± 0.15 MeV KEDR 
aims to obtain a final accuracy of 0. 15 MeV, and BESlll aims at a precision better than 0. 1 MeV 

2.2 Tau lifetime 

The tau lepton lifetime is known up to 0.3% and is the least precisely known factor in the 



couphng ratio expressions ^ and IZ4 BaBmi. has presented a preliminary measurement of the tau 
hfetime with an error comparable to the present world average: 

= 289.4 ± 0.9 ± 0.9 fs [0]. (2.6) 

The measurement uses about 80 fb^' of data and is based on an extremely pure (99.4%) yet scarcely 
efficient (0.2%) selection of 1 against 3-prong events in the cm. system, where the 1-prong track is 
an identified electron. Electron identification is used because it is more efficient and less contam- 
inated with hadrons with respect to muon tagging. The selected tau candidates are about 300,000 
and include about 0.2% hadronic background, 0.4% Bhabha background and a negligible amount 
of two-photon events. 

The measurement is based on the reconstruction of the decay length of the tau that decayed 
into the 3-prong tracks. Using a novel technique aimed at minimizing the systematic dependence 
on the detector alignment, the tau decay vertex is first computed in a plane transverse with respect 
to the beam axis. The transverse decay length is computed within the transverse plane by projecting 
the vector from the luminous region center to the tau decay vertex along 3-prong total momentum 
direction (which approximates the tau flight direction). The tau decay length is finally reconstructed 
by projecting the transverse decay length onto the 3-prong total momentum direction. 

The mean decay length is determined with an average, abstaining on purpose from weighting 
events according to their estimated errors in order to minimize systematic effects from alignment 
and detector material modeling. The mean lifetime is determined using the Monte Carlo prediction 
of the average tau momentum, using the KKMC generator [[l^], which includes complete 2nd order 



radiative corrections. The measurement offset that originates from tracking errors correlations [ |1 1[ ] 
and from approximating the tau momentum direction with that of the 3-prong total momentum is 
subtracted using the Monte Carlo simulation. 

Finally, the contribution of background is subtracted. While hadronic background is simulated, 
a statistically adequate Monte Carlo simulation of Bhabha events is impractical because the relevant 
cross-section is about 20 times the tau production cross section (ss 1 nb), therefore a data control 
sample is used to estimate both the Bhabha contamination and its decay length distribution. 

Systematic uncertainties come mainly from the reliability of the measurement bias subtraction 
using Monte Carlo, from detector alignment, from the mean tau momentum Monte Carlo sim- 
ulation, and from background subtraction. This measurement includes a study of the effects of 
detector misalignment. Decay length shifts with respect to a perfectly aligned detector are mea- 
sured on simulated Monte Carlo events by refitting tracks from coordinates taken on a detector that 
is purposefully distorted. Six representative distortions are applied by displacing the silicon vertex 
detector wafer positions according to the observed distortions and uncertainties that are derived 
from reconstructed data. 
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Figure 2: Selected tau lifetime measurements (left) and check of the Standard model (SM) prediction of 
universal leptonic couplings to the W (right) combining the present tau lifetime world average with the 
BaRar preliminary result. The thickness of the oblique line represents the uncertainty of the SM constraint, 
and is dominated by the uncertainty on the tau mass. 



This preliminary measurement contributes the largest experimental improvement in recent 



years for the coupling ratio expressions 2.1 and 2.2. Figure H reports the measurement compared 



with the present world average and previous selected measurements ||13|]. Combining the BaBj^ 



2004 result with the present world average assuming no systematic error correlations we obtain 

= 290.15 ±0.79 fs. (2.7) 

Using the present world averages, we present an updated check of lepton universality in Figure ^ 
and updated determinations of the coupling ratios 

^ =0.9982 ±0.0020, —= 0.9980 ±0.0020, ^ = 0.9981 ±0.0017, (2.8) 

gt gt gt 

where ge^^ is determined assuming = g^ holds for the theory. 

In the recent past, LEP experiments improved considerably the experimental precision on the 
tau lifetime profiting from ideal conditions in most respects but statistics: there high momentum 
tracks had small impact parameter errors due to multiple scattering, tau events had a distinctive 
topology that permitted a pure and efficient selection against backgrounds, vertex detectors pro- 
vided precise tracking close to the origin and systematic uncertainties from detector misalignment 
were reduced thanks to the complete and uniform acceptance in the azimuthal angle [[Til]. B- 
factories appear to be the only facilities where the tau lifetime measurement can be improved in the 
near future and they can overcome with statistics the disadvantages related to increased multiple 
scattering due to lower momenta and to less favourable physics conditions for an efficient and pure 
selection. 

2.3 Tau leptonic branching fractions 

B-factories have not yet succeeded in matching LEP experiments regarding the measurement 
of the tau leptonic branching fractions, because the systematic precision is limited by uncertainties 
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on the normalization of the number of produced tau leptons, arising from uncertainties on the 
integrated luminosity corresponding to the analyzed event samples and on the cross-section for tau 
pair production at the T(45'). Uncertainties on luminosity are about 1% to be compared with better 
than 0.1% for LEP experiments. Until recently, the uncertainty on the tau pair cross-section was 



estimated by comparing the KoralB [14] and KKMC predictions to be 2.2%. However there has 



been remarkable progress that led to a new estimate of the tau pair cross section at the T(45) : 

a(f'+e" ^T+T") =0.919 ±0.003 nb [|]. (2.9) 



Furthermore, B-factories may measure ratios of branching fractions (see eq. [2.3[ ) whose uncer- 
tainties won't be limited by the understanding of luminosity and cross-section, but only by the 
systematics related to electron versus muon identification and to background suppression and sub- 
traction. 



3. Lepton Flavour Violation tests at the B-factories 

3.1 Common analysis features 

The typical tau LFV decay search at the B-factories selects low track multiplicity events that 
have 1 against 1 or 3 tracks in the cm. frame. The thrust axis is used to define two hemispheres, 
each of which is then examined for consistency with a tau LVF decay, while the other one must be 
compatible with a known tau decay. Unlike known tau decays, which include at least one neutrino, 
the reconstructed products of a LFV tau decay are expected to match the tau mass and half the 
cm. energy within the experimental resolution. It is worth noting that physics effects also limit the 
experimental accuracy in reconstructing the parent tau energy and mass from its decay products: 
initial and final state radiation affect the tau energy itself before decay, and radiation in decay 
and Bremsstrahlung from the decay products change the reconstructed energy and invariant mass. 
The energy is reconstructed with a typical resolution of 50MeV and, when using a total energy 
constraint to half the cm. energy, the invariant mass is reconstructed with a resolution of about 
lOMeV. Selected events around the expected energy and mass within 2 or 3 standard deviations 
are then investigated looking for an excess over the expected background. 

The amount of expected background is normally estimated using the distribution shapes from 
the Monte Carlo simulation normalized to the observed events in a two-dimensional sideband re- 
gion around the signal region on the energy-mass plane. The signal efficiency is estimated with a 
Monte Carlo simulation and typically lies between 2% and 10% depending on the channel. Typical 
cumulative efficiency components include 90% for trigger, 70% for geometrical acceptance and 
reconstruction in the detector, 70% for reconstructing the selected track topology, 50% for particle 
identification, 50% for additional selection requirements before checking the reconstructed energy 
and mass, and 50% for requiring consistency with the expected energy and mass. The selection effi- 
ciency and background suppression are optimized to give the best "expected upper limit" assuming 
that the data contain no LFV signal. The optimization and all systematic studies are completed 
while maintaining the experimenter "blind" to data events in the signal box in the energy-mass 
plane, in order to avoid experimenter biases. 

When the expected background in the signal region is of order one or less, the number of signal 
events is normally set to the number of observed events minus the background, while in presence 
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of sizable background the numbers of background and signal events are concurrently determined 
from a fit to the mass distribution of events that have total energy compatible with the expected 
one. 



3.2 Results 

LFV decays can be grouped in the following categories: tau to lepton-photon (T £7, where 
I = e,}JL), tau to three leptons or one lepton and two charged hadrons (t ^i^2^3> I' ihihz), tau 
to a lepton and a neutral hadron (t — > ih'^', where h'^' = 7i^,ri,ri',K^,etc.). 



Belle has reported preliminary results for T — > /xy and z ej [15] based on the analysis of 
535 fb~^ of data. The T jxy (t — > ey) searches have a 5.1% (3%) signal efficiency within a 2a 
elliptical signal region in the energy-mass plane plane. A two-dimensional unbinned extended max- 
imum likelihood fit for signal and background in the signal region obtains -3.9l^ j (-0-14+2:45) 
signal and 13.9+4 ^ (5.14+2 si) background events. Frequentist 90% CL limits are obtained by run- 
ning Monte Carlo simulations in which the signal is increased until 90% of the fits obtain more 
than the observed signal events: =^(t piy) < 0.45 • 10"^ and =^(t ey) < 1.2 • 10"'. 



BaBar has recently submitted for publication [ ]16| ] and Belle has recently reported [ |17| ] im- 
proved results on tau to three leptons LFV searches based on enlarged event samples of 400 fb"^ 
and 535 fb"^ respectively. The two analyses have similar signal efficiencies (5.5%— 12.5%) al- 
though different strategies are adopted to define the signal regions in the energy-mass plane: Belle 
uses ellipses large enough to contain 90% of signal events, while BaBar uses rectangular signal 
boxes optimized to obtain the lowest expected upper limit in case there is no signal. The Belle 
selection is significantly more effective in suppressing the background, which is expected to be 
0.01—0.07 events in all channels but the three electron one, where 0.4 events are expected because 
of Bhabha contamination. The background in the signal region is estimated from the mass distri- 
bution sidebands, assuming it is constant, using looser selection criteria to get reasonable samples 
close to the signal region. No details are given on how the extrapolation is done for the full selec- 
tion. BaBj^ expects 0.3—1.3 background events, i.e. of order one, since the selection is optimized 
for the best expected upper limit, at the risk of getting background events in the signal region. 
Belle observes no candidate signal events in 535 fb"' of data in all modes, and calculates upper 
limits using Feldman and Cousin ordering 18] in the range [2.0—4.1] • 10"^, depending on the 



mode. BaBar observes from to 2 events in 376 fb of data, and calculates upper limits accord- 



ing to Cousin and Highland prescription [19| with no Feldman and Cousin ordering in the range 
[4-8]-10-^ 

BaBar has recently published new results on tau LFV decays into a lepton and a hadron pseu- 



doscalar tt",!],!]' [20|. In these analyses both of the rj ^ yy and the T] ^ 3 TT decay modes are 
used for the analyses, and 77' candidates decaying both to rjln and yin are considered. The ex- 
pected background per channel is between 0.1 and 0.3 events. Summing over all ten modes, 3.1 
background events are expected, and 2 events are observed. 

Belle has recently reported improved results on tau LFV decays into a lepton and a vector 
meson V° pill, with = (j), CO, K**^ or K*^\ using 543 fb^^ of data. No excess of signal events 
over the expected background is observed, and upper limits the branching fractions are obtained in 
the range (0.7-1.8) • lO"'' at 90% CL. 
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Channel 


Belle 




BABm 


combined 


UL90 


Lumi 


UL90 


Lumi 


UL90 


Lumi 




(10-^) 


(fb-i) 








\W ) 


M7 


0.5* 


535 


7 


232 


16 


1 yj 1 


ey 


1.2* 


535 


1 . 1 






if\i 




0.65* 


401 




339 


51 


140 


}xr]' 


1.3* 


401 


1 3 


339 


53 


740 


er\ 


0.92* 


401 


1.6 


339 


0.45 


740 


erj' 


1.6* 


401 


2.4 


339 


90 


740 




1.2* 


401 




339 


58 


740 




0.8* 


401 


1 3 

J- 


339 


44 


740 


m 


0.20-0.41* 


535 


0.4-0.8* 


376 






£hh' 


2-16 


158 


1-5 


221 






^\/° 


0.7-1.8* 


543 










^lKs 


0.49 


281 










eKs 


0.56 


281 










An, An 






5.8-5.9* 


237 






AK,AK 






7.2-15* 


237 












40-89 


211 







(* preliminary) 



Table 1: Summary of 90% CL upper limits on tau LFV decays from the B-factories. An asterisk indicates 
a preliminary result, h and h' denote a charged pion or kaon. Banerjee's combination of a subset of these 
channels is also included. 



BaB^ reported also on less conventional searches of tau LFV decays into A7r [26| and of LFV 
in tau production {e^e^ It) [27], finding no signal. 



The above results and additional B-factories LFV results [^-^] are summarized in Table 
At the Tau2006 conference in Pisa, Swagato Banerjee presented frequentist combinations [ p5| ] of 
some measurements, which are also included in the table. 

3.3 Prospects 

While Belle plans to run until it will collect 1 ab^^ of data, BAB/R is funded to run until 
September 2008, when it expects to reach a total integrated luminosity of about 0.8 ab^^. In case 
there is no signal, the expected upper limits on the number of selected signal events will improve 
depending of the amount of irreducible background in each channel: 

• when the expected background is large (A'^bkg ^ 1). the expected upper limit is N^^ ^ 
1.64VA^bkg; 

• when the expected background is small (A^bkg ^ 1)> using [g one gets ^ 2.4. 



Reducing the background below few events does not much improve the expected limit if significant 
efficiency is lost in the process, therefore optimized searches often enlarge the acceptance until 
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A^BKG ~ 1- For the cleaner channels, analyses can be optimized for an increased data sample to 
keep A^BKG ~ 1 without loosing a significant part of the signal efficiency: in this best case scenario, 
the expected upper limits will scale as A^bkg/-^^ i-C- as 1 /^. On the other hand, if no optimization 
is possible, just keeping the current analyses will provide upper limits that scale as a/A^bkg/-^> i-C- 
as \/^/^. 

Mike Roney has recently estimated [ ^8| ] that for T ^ ly analyses the background coming from 
T Ivv decays and initial and final state photons can be considered "irreducible" at its present 
level (20% of the total). In these conditions, the final combined Belle and BfiBPR data set will allow 
for expected upper limits in the range [0.1—0.2] • 10^^. The actual combined upper limit obtained 
by Banerjee [ p^ ] for T ^ /^y is already in that range as a consequence of a downward fluctuation 
of the observed events with respect to the expected background. 

Other tau LFV decay channels (t Ul, T — > ihh, z — > £h^) do not yet appear to be background 
limited and their expected expected upper limits with the full B-factories dataset are about 0.1- 
10-^. 

Beyond the current B-factories facilities, there are proposals for super B-factories [29] that 
would permit a 100 fold increase in the size of the tau pairs sample: this would allow probing tau 
LFV decays at the lO^'^-lO^i^ level. 
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